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ABSTRACT 

We present near infrared data on the luminous interacting system Mkn 266 
(NGC5256), comprising 2 /mi continuum, and Br7 and 1-0 S(l) emission line images, 
together with K-band spectra. We have fitted stellar templates to the continuum, 
allowing us to account for all the stellar features and hence detect even faint gas 
excitation emission lines, including 8 and 11 H2 lines in the SW and NE nuclei 
respectively. Population diagrams for the excited H2 molecules indicate that most of 
the 1-0 S(l) in each of the nuclei has a thermal origin. We discuss this with reference 
to the observed morphologies, especially that of the 1-0 S(l) line. In particular, the 
core of 1-0 S(l) in the NE nucleus is more compact than the 2 /jm continuum; while in 
the SW nucleus the 1-0 S(l) is significantly offset by 500 pc from the continuum (and 
other) emission. Lastly we address the issue of the region midway between the two 
nuclei, where previously a strong source of radio continuum has been observed. These 
results are set in the context of interacting galaxies where shock excited emission might 
be expected to occur as a direct consequence of the interaction. 

Subject headings: galaxies: individual (Mkn 266) - galaxies: interactions - galaxies: 
nuclei - infrared: galaxies - line: formation 
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1. Introduction 

Prompted by the discovery from CO observations that the bulk of the molecular gas in 
close mergers can quickly congregate between the nuclei (Gao et al. 1997, Scoville et al. 1997), 
we have begun to look at number of such systems, imaging the 1-0 S(l) and Br7 lines since 
the morphologies of these with respect to the continuum and each other can provide important 
information. The rationale is that while star formation models allow a quantitative analysis of the 
observed fluxes, the results can be quite misleading if the distribution of the emitting regions is not 
considered. Compare, for example, an object in which the Br7 and 2 fim continuum come from the 
same region (in which a single burst of star formation is possible) , to one in which they are offset 
from each other, requiring two distinct epochs of activity. A classic example is given by Gonzalez 
Delgado et al. (1997) who found both Wolf-Rayet and Can features in the super star cluster A of 
NGC 1569, implying the simultaneous presence of hot massive stars and red supergiants. Their 
intepretation required the existence of non-coeval stellar populations in the same stellar cluster. 
The difficulty was resolved when De Marchi et al. (1997) found from HST imaging that this was 
actually 2 distinct clusters. Such high resolution is not always required, and in a small sample 
of blue compact dwarf galaxies observed by Davies et al. (1998), the Br7, 1-0 S(l), and 2/xm 
continuum all had strikingly different morphologies on scales of a few arcsec (~100pc). In such 
cases a simple interpretation of the emission in terms of just star formation is not valid. 

Another issue which has posed a significant barrier to extragalactic work is the detection of 
the fainter H2 lines, which is essential to determine the fractions of H2 which are thermally and 
non-thermally excited. Previously, except in a few cases such as NGC 6240 (Sugai et al. 1997b) 
which has extremely bright 1-0 S(l), this has typically relied on very few line ratios, often given 
only as upper limits and extracted from rather large apertures (eg Goldader et al. 1997). It is 
crucial to detect lines from the v=2 (or higher) bands: a single thermal model can nearly always 
provide a reasonable fit to lines from the v=l-0 transitions. It is also important to have more 
than one v=2-l line, since ratios involving only a single such transition can be ambiguous: for 
temperatures in the range 1500-2500 K, the thermal 2-1 S(l)/l-0 S(l) ratio varies from 0.03 to 
0.15. Thus if a ratio (or upper limit) of 0.15 is observed, while the 1-0 S(l) is constrained to be 
> 75% thermal, the total H2 cooling can only be attributed to somewhere in the range 35-100% 
thermal depending solely on a fairly small change in thermal excitation temperature. With the 
increasing sensitivity of near-infrared detectors it is now possible to use small apertures and obtain 
the high signal-to-noise spectra required to detect higher vibrational lines. Additionally, fitting 
the continuum with stellar templates, rather than using a power-law or polynomial fit, is necessary 
in order to take account of the stellar features and allow weak lines from the excited gas, such as 
H2 lines, to be measured reliably. 

We address these issues in this paper, in which we present a detailed analysis of near 
infrared emission line images and spectra. We focus on Mkn266, a luminous infrared galaxy 
(LrR = 3 x 10 11 L Q ) at a distance of 115 Mpc (1" = 560 pc). It has 2 prominent nuclei 10" apart, 
with a common envelope, suggestive of a merger (Mazarella et al. 1988, hereafter MGAH). There 
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have been very few detailed studies of Mkn 266, the two most decisive being those of MGAH and 
Wang et al (1997). The former authors presented an Ha map which revealed arcs that may be 
tidal tails or remnant spiral structure. They also mapped the radio continuum, which appears 
similar to the 2 ^m continuum but with a strong region of extended emission halfway between the 
nuclei. The latter study found evidence for a superwind in their X-ray imaging and spectroscopic 
data. 



2. Observations and Data Reduction 

2.1. Fabry-Perot Images 

Images of the 2 ^m continuum, and 1-0 S(l) and Br7 emission lines were obtained with 
the UKIRT 3.8-m telescope on Mauna Kea, on the nights of 12 and 13 May 1998. A K-band 
Fabry-Perot etalon (FP), set in a collimated beam, was used in conjunction with IRCAM3 (a 
256 x 256 InSb array) at the Cassegrain focus. The scale of 0.286" per pixel gave an unvignetted 
field of view of ~60". Order sorting was achieved with cooled narrow band filters: 2.4% FWHM 
for the Br7 line, and 4.25% FWHM for the 1-0 S(l) line. The latter filter is normally used for 
continuum imaging, but was required due to the large recession velocity (> 8000 km s -1 ) of the 
galaxy which shifted the line out of the standard filter pass-bands. The 325 km s" 1 spectral 
resolution of the FP and the wavelength calibration were checked several times during the night 
using the 2.1171 /xm line from a Krypton lamp, permitting adjustments for temperature changes 
(approx. 40kms~ 1 /°C). 

Due to the relatively small size of Mkn 266 we imaged it at 2 positions on the detector, 
symmetrically off-axis from the centre of the FP (14" North- West and 14" South-East) . The 
transmitted wavelength varies off-axis as S\ oc r 2 , reaching SX = 160 km s -1 at r = 30", hence 
the difference is only 5X = 40kms~ 1 at r = 14". This was taken into account when setting the 
FP wavelength. Other concerns include the small 15kms _1 offset for the velocity of the earth 
relative to the sun in the direction of the object on that date, which was not included, and the 
uncertainty in the heliocentric line velocity. Our adopted value was 8385kms _1 , the average for 
the two nuclei determined from optical spectra by Veilleux et al. (1995), resulting in additional 
offsets of ±25kms~ 1 for each nucleus. Although each of these adjustments is small, they must 
all be considered: with a FWHM resolution for the FP of 325 km s -1 , a total offset of 100 km s -1 
will mean that the line flux is underestimated by 35%, but a reduction to only 50kms _1 (similar 
to the maximum we expect) results in a net effect on the line flux of less than 10%. A similar 
consideration applies to the line dispersion, which can reduce the line flux detected if it is much 
greater than 100 km s -1 . Typical corrections are 20% for 200kms _1 , and 35% for 300kms _1 , but 
it should be noted that they combine with velocity offsets non-linearly. Line dispersion effects 
were not included in the line flux estimation since the line-widths are not known, although they 
are expected to be relatively narrow for starbursts. 



-4 - 



Observations were carried out in blocks, each of which had 4 integrations of 3 minutes per 
pointing, consisting of: 2 at the on-line wavelength, and one each for continuum at slightly 
shorter and longer wavelengths (shifted with respect to the line centre by -800 and +1000 km s -1 
respectively). Standard A-type stars (BS4344 and BS5142) were observed between each set of 2 
or 3 blocks. The total on-line integrations were 48mins for Br7, and 120mins for 1-0 S(l). 

IRAF was used for all image processing and analysis as described below. Flatfields were made 
using opposite positions after dark current subtraction. These consist of the camera flatfield, and 
a smooth circularly symmetric contribution from the FP. Due to small temperature variations, the 
latter changed slightly between successive integrations (at the 1% level), leaving a residual ripple. 
Although this makes little difference to the continuum counts, it is crucial to the detection of line 
flux since the peak line counts per pixel are only 2% of the sky background. Thus after each frame 
was dark subtracted and flatfielded, it was also divided by the median filtered combination of itself 
rotated every 22.5° about the centre of the ripple pattern. A constant value for the sky could then 
be subtracted (rather than using sky frames which would increase the noise). The frames were 
aligned on the NE nucleus since it is more compact than the SW one, but this could introduce 
systematic offsets due to the line emission in the on-line frames. In order to check that there was 
no such bias, the centroid positions between pairs of successive off-line and on-line images (at the 
same telescope pointing) were compared. For the Br7 line these were x = — 0.17 ± 0.47 pixels and 
y = 0.05 ± 0.27 pixels; for 1-0 S(l) they were x = -0.10 ± 0.12 pixels and y = 0.24 ± 0.20 pixels 
(1 pixel = 0.286"). For both lines, the mean difference was less than the RMS variation and 
systematic effects should not be more than ~0.l". 

Flux calibration was achieved with the standard stars BS4344 (type A4V, K=6.30) and 
BS5142 (type A3V, K=5.18). The frames were treated in the same way except that sky frames 
were subtracted and there was no need to correct for the ripple pattern. The profiles of the 
order sorting filters were included, making as much as 15% difference from assuming a simple 
boxcar shape due to secondary FP transmission orders in the wings of the pass-band. The image 
resolution was determined from the stars to be 0.77" for the 1-0 S(l) line and 1.34" for the Br7 
line. 

The continuum image was deconvolved with 30 iterations of the Lucy algorithm implemented 
in IRAF, after halving the pixel size. The iterations were stopped once false structure due to noise 
peaks began to show in the fainter emission. At this point the stronger emission peaks show only 
genuine features; for example, the asymmetry in the NE nucleus in the deconvolved image is also 
apparent (just) in the direct image. The resolution of 0.36" was determined by including at one 
edge of the image an extra region with the same noise and a number of Gaussians convolved with 
the PSF; after deconvolution, the recovered and original sizes were compared. 
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2.2. Longslit Spectroscopy 

Spectra of Mkn 266 were obtained as part of the UKIRT service observing programme on May 
1 1999 using the CGS4 spectrometer; the orientation of the slit, 33° East of North, is indicated on 
Fig. |l]. The total integration time was 64 min. Frames are sky subtracted, flatfielded, and co-added 
on-line using CGS4DR; further data reduction was carried out with IRAF. The dispersion axis was 
aligned to pixel rows, and the spatial axis to pixel columns with reference to an argon arclamp and 
hence simultaneously calibrating the wavelength scale. Standard A stars (BS 5023 and HD 105601) 
were used for flux calibration and atmospheric correction after interpolating over Br7 and Paa 
with Lorentzian profiles and making a small correction to the same airmass as Mkn 266. The 
spatial resolution was determined from the standard stars to be 1.06", well matched to 2-pixel 
sampling. Integrated spectra of each nucleus were extracted in order to look at weak features, and 
additionally the spatial extent of the brighter features was considered. The heliocentric velocities 
of the two nuclei were measured as 8440 km s" 1 (NE) and 8390 km s" 1 (SW), in good agreement 
with those of Veilleux et al. (1995). 

3. Emission Line Images 

The continuum image (upper left and, deconvolved, upper right) in Fig [l] shows that the 2 
nuclei are separated by 10.5", 5.9 kpc at the distance of the galaxy. The NE nucleus is resolved 
with a size of 0.8" xO. 6". As Fig |2| shows the 1-0 S(l) is more compact than this: comparison 
via quadrature correction with the continuum would put its size at <J 0.45". This is similar to 
the extent of the radio source observed at 6 cm by MGAH, which was marginally resolved at a 
resolution of 0.3" x 0.4", and found to have a spectral index of 0.61±0.04. 

The continuum in the SW nucleus has a morphology very similar to the radio continuum: 
central source with extensions to ~1" at position angles slightly west of North and east of 
South. The Br7 may also include these regions but the poorer resolution rules out a definitive 
statement. Nevertheless, the presence of a fairly substantive continuum component argues against 
the suggestion of MGAH that these are radio jets from the Seyfert nucleus, and tends towards 
interpretation in terms of star formation. An interesting point is that the high CO index measured 
in the spectrum by Goldader et al. (1997), CO p h = 0.16 ± 0.03, suggests the presence of late-type 
supergiants, which would only be expected in a young stellar population. Such a conclusion is 
consistent with our spectra which suggest a 50% contribution from supergiants in the K-band, 
twice that of the NE nucleus. The implied spatial scales on the order of 1 kpc are consistent with 
what is known about circumnuclear star formation occuring at inner Lindblad resonances around 
an AGN/starburst. It should be noted, however, that there is currently no evidence for the barred 
potential in a disk which would lead to gas inflow to such resonances. Indeed, the Ha and optical 
morphologies (MGAH) are highly perturbed. 

The 1-0 S(l) in this region emission is rather unexpected, the peak being offset from the 
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continuum by 0.9". As described in Section although the images were aligned on the NE 
nucleus, systematic effects cannot account for these offsets. The 1-0 S(l) line map has a higher 
signal-to-noise and good resolution, so we can be confident that the offset is real. For the Br7 map 
a big uncertainty lies in the low signal-to-noise (peak pixel in the smoothed image is 7.5a), and 
the relatively poor resolution (1.3"). We therefore cannot claim any significance in the Br7 offset. 
Astrometry with the 6 cm radio map of MGAH suggests that the Br7 and 2 /jm continuum are 
aligned with the Seyfert nucleus. 

4. Spectra 

An integrated spectrum of each nucleus was extracted, covering 5 pixels (3.05") along the slit 
by the slit width of 1.22", and these are shown in Fig. ||[ Additionally, the spatial extent of some 
stronger features could be mapped, as shown in Fig. ||. Clearly the system is far from trivial, and 
it would be useful to analyse all the line ratios with high spatial resolution, something which may 
be possible with adaptive optics and integral field spectroscopy in the (near) future (eg see Davies 
et al. 1999). Here we present an analysis of the detailed integrated spectra together with some 
spatial information from both spectra and images. 

4.1. Continuum Fitting & CO Absorption 

In order to recover emission lines which are either faint or superimposed on absorption 
features it is important to subtract an accurate continuum rather than fit a simple line or curve 
as is often done. This has the double benefit of also providing a qualitative estimate of the stellar 
population contributing to the 2/im emission. We used a minimum chi-square fitting routine with 
the following restrictions: 

(1) only regions of the continuum distinct from emission lines were included in the chi-square 
estimations; 

(2) a single value for extinction was fitted; 

(3) a hot dust component (eg for an AGN) was not included as it did not improve the fits. 

The CO absorption indicates that late-type giant & supergiant stars dominate the continuum 
at these wavelengths and so we used generic spectra of these types, created by adding several 
examples of each from the library of Forster-Schreiber (1999) and convolving them to the same 
resolution as our data. The individual templates have varying spectral coverages (dependent on 
the resolution), resulting in a trade-off between large coverage of 2.00-2.44/iin or a wider variety 
of stellar types but covering only 2.27-2.40/iin. The latter range of templates includes MI, KI, 
KOIII, K5III, MO III, which cover the range of types required; adding any other available spectra 
(eg between K0 and K5 giants) would not provide any further useful information about the stellar 
population. The former range of types misses out the MO III and, more seriously, KI. Nevertheless, 
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both approaches gave the same qualitative result, and so we used the former one in order to 
provide the longest possible baseline for subtraction. The stellar types used were MI, KOIII, and 
K5III, and the best fits had RMS errors of ~0.06 per data-point (for the scales in Fig. |3|). The 
baseline of the fitted continuum was extended down to 1.85/zm by adding a much lower resolution 
spectrum of a late-type giant star, also reddened with the same extinction, from the PEGASE 
database (Fioc & Rocca-Volmerange 1997). 

The result showed that about 20-25% of the 2/im continuum in the NE nucleus is from 
late-type supergiants, while the rest is from giants, and the extinction is Ay ~ 3.2 mag; in the SW 
nucleus about 50% is from supergiants, with Ay ~ 4.6 mag. 

We have additionally tried to constrain the stellar population by measuring the CO index 
directly. We have used two indices, CO' of Goldader et al. (1997) in the 2.30-2.34 /mi interval, 
and CO sp of Doyon et al. (1994) in the 2.31-2.40 /im range. Converting these to CO p h, using 
the relations given by Goldader et al., gives similar values for the two indices. These values are 
COph = 0.12 (NE) and 0.16 (SW). The latter value verifies that there is a higher supergiant 
fraction in the SW nucleus, and is the same as that measured by Goldader et al. even though they 
used a much larger 3"x9" aperture. 

4.2. H2 Excitation 

The process described above removes much of the structure in the continuum, including such 
features as Mg (2.14/xm), Na (2.21/xm), Ca (2.26/im), etc. It allows not only weak lines to be 
measured, but also the Q-branch H2 lines - which is useful since the 1-0 Q(3) and 1-0 S(l) lines 
come from the same upper level, as do 1-0 Q(2) and 1-0 S(0), providing a consistency check when 
plotting a population diagram for the molecules. 

The H2 line ratios relative to 1-0 S(l) in each nucleus are given in Table ||, and level 
population diagrams are shown in Figs. || and ||, which have been derived assuming the local 
thermal equilibrium ortho-para ratio of 3. This should be valid for thermally excited H2, but it 
is generally found that for H2 excited by UV fluorescence the ortho-para ratio is ~2. Sternberg 
& Neufeld (1999) argue that this is because the ortho absorption lines become optically thick 
sooner than the para absorption lines, so UV pumping of the ortho states is less efficient and 
the ortho-para ratio appears to be reduced to a limiting case 

of ~ -v/3 = 1.7. That is, the UV 
excited states (those observed) have an apparent ortho-para ratio of 1.7, while the total H2 
ratio remains at 3. This effect only occurs if the absorption lines become optically thick to UV 
pumping, and is included in the line strength calculations of Black & van Dishoeck (1987) from 
whom we take our models of UV fluorescence. However, for the purposes of fitting fluorescent 
models to the population diagram, the ortho-para ratio used does not actually matter as long 
as the same ratio is used for both the observed and modelled spectra when converting line 
strengths to level populations. For simplicity, therefore, since we do not know a priori the relative 
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thermal/non-thermal contributions (which is what determines the observed ortho-para ratio), we 
have adopted the same ratio as in the LTE case. 

For one model, we have fitted a combination of thermal and non-thermal components, 
with the 3 free parameters being the relative non-thermal contribution /uv to the 1-0 S(l) line, 
the temperature T of the thermal component, and the absolute scaling. For the non-thermal 
component we used fluorescent model 14 from Black & van Dishoeck (1987), although this 
could equally well be formation pumping which results in a similar spectrum. For the thermal 
component we used the Boltzmann distribution at temperature T given by 

^ oc e- E ^ kT with N u oc — 
9 Ail 

where N u is the column density of H2 in the upper level of excitation energy E n , of an observed line 
which has wavelength A, flux F, and transition probability Ail (quantities are listed in Table ||); 
g is the degeneracy given by the product of the rotational and spin degeneracies gj x g$ where 
gj = 2 J + 1 and gs = 3 for odd J and 1 for even J. In the optically thin regime, these equations 
can be used to estimate the mass of hot gas Mh 2 from the line luminosity L since 

L> — Ju All "T 

A toh 2 

where / u is the fraction of hot H2 in the upper level of the line , and mn 2 is the mass of a hydrogen 
molecule. Using the normalisation from Scoville et al. (1982) for gas thermalised at 2000 K that 
for the v = 1, J = 3 level / u = 0.0122, the hot gas mass can be estimated from the 1-0 S(l) line 
luminosity as 

£i-os(i)(£©) = 62OM H2 (M ) 

In the population diagram, parameters were derived by minimising chi-square; confidence 
regions were determined by systematically varying each of T and /uv (which are not necessarily 
independent) while simultaneously optimising the absolute scaling, until chi-square had increased 
by one standard deviation. For the NE nucleus T = 1500 ± 60 K and /uv = 0.19 ± 0.02; for the 
SW nucleus T = 2460 ± 410 K and / UV = 0.29 ± 0.12. 

For a second model we attempted to fit two thermal components, since it is not implicitly 
ruled out by the data. For the NE nucleus, the reduced chi-square xt = 1-96, rather larger than 
the expectation value of (0.91) for 7 degrees of freedom. As well as this statistical arguement, there 
is a physical reason against the model: the temperatures derived are 1200 K and 5260 K. While 
the former only varies by ±100 K, the la limits on the latter encompass the range 4000-7600 K. 
Since H2 rapidly dissociates at temperatures above 4000 K, it is not possible to get a strong 
thermal spectrum at the temperatures indicated here. Thus we can rule out this model for the 
NE nucleus. For the SW nucleus, the two-component thermal model effectively reduces to a 
single thermal component which is itself simply a subset of the previous thermal plus non-thermal 
model. The value xi = 1-95 for a single thermal component supports the result of the thermal 
plus non-thermal model above that a pure thermal model is unlikely. 
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4-2.1. UV Fluorescence 

If the non-thermal H2 component arises in photo-dominated regions due to UV fluorescence at 
the edge of Hn regions, the models of Puxley et al (1990) can be used to constrain the compactness 
of the star forming regions. The diagnostic needed is the fluorescent 1-0S(1)/Br7 ratio, which we 
estimate to be 0.35 ± 0.04 (NE) and 0.16 ± 0.07 (SW); the latter ratio may be intrinsically higher 
if any of the Br7 is related to the AGN. 

The ratios derived are rather higher than expected, and the models suggest that they are 
typical of star formation which is either fairly diffuse or relatively evolved (so the most massive 
ionising stars are no longer present), both of which seem rather unlikely. The former option 
contradicts observations that intense star formation appears to occur mostly in compact clusters 
2-3 pc across with masses up to 10 6 M Q (eg Meurer et al. 1995) which would have ionising fluxes 
on the order of 10 52 ~10 53 sec -1 and hence, according to the model, rather small ratios. Indeed, 
Davies et al. (1998) found that the 1-0S(1)/Br7 ratios of such clusters, averaged over the central 
few tens of parsecs in the nuclei of blue compact dwarf galaxies, were <^0.1. The latter option 
can be refuted by considering the Br7 equivalent width, Wbi^, which implies a maximum age of 
50 Myr for the clusters dominating the 2 fim continuum. Since this time span is only greater than 
the main sequence lifetime of stars more massive than ~15M , many highly ionising stars will 
still be present, implying a low UV flourescent 1-0S(1)/Br7 ratio. The high Hei/Br7 ratios of 
0.62 ±0.04 and 0.47 ±0.03 in the NE and SW nuclei respectively, also point towards very hot stars. 
Although not a unique indicator of stellar temperature (due to resonant and, at high densities, 
collisional effects, Shields 1993), such high values suggest T e f f = 30-40 xl0 3 K equivalent to stars 
of mass 15-30 M Q . As discussed below, some other process as well as UV fluorescence is probably 
needed to account for much of the non-thermal H2 emission. 

At high densities, ^> 10 5 cm -3 , collisional effects will tend to thermalise the H2 molecules, 
altering the emission spectrum. This means that in principle UV fluorescence can also give rise to 
a spectrum that is apparently thermal, an effect used to explain the anti-correlation between the 
ratio of 2-1 S(l)/l-0 S(l) and the intensity of 1-0 S(l) (Usuda et al. 1996). Here, since it appears 
already that fluorescence can account for only part of the weak non-thermal 1-0 S(l) emission, it 
is unlikely to contribute more than a few percent to the much stronger thermal part. 

4-2.2. X-ray Excitation 

In order to determine whether X-ray irradiation can be a potential source of H2 line excitation 
we need to estimate the 1-100 keV X-ray flux. The 0.1-2.0 keV X-ray luminosity has been 
determined by Wang et al. (1997) who analysed the X-ray emission from Mkn266, making use 
of both the High Resolution Imager (HRI) and Position Sensitive Proportional Counter (PSPC) 
instruments on RO SAT to identify a soft diffuse halo as well as 3 harder compact sources. They 
found that one associated with the SW nucleus has a luminosity of 1.3-7.5 x 10 41 ergss _1 and 
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the NE nucleus has 0.9-5.1 x 10 ergss , depending on the absorbing hydrogen column density 
(iVn = 0.016 &: 0.3 x 10 22 cm -2 ). However, they found that only thermal plasma models fit the 
data satisfactorily, an unusual result given that the SW nucleus hosts an AGN for which the 
canonical hard X-ray spectrum is a power-law F v oc v a with index a = —0.7. That there is no way 
to be sure which sources seen with the HRI correspond to which parts of the spectrum observed 
with the PSPC adds confusion. Since the SW nucleus contributes <^ 20% of the total counts, these 
may not have much impact on the 0.1-2 keV spectrum, and it is possible that the spectrum at 
energies lower than 2keV is dominated by star-formation, while a Seyfert-like power-law would 
only emerge at higher energies. This is seen in some Seyfert2s, such as Mkn3 (Serlemitsos, Ptak 
& Yaqoob 1996). 

This could well the case for the SW nucleus, where we make the reasonable assumption that 
the AGN dominates the X-ray luminosity, so we can put very comfortable upper limits on the 
extrapolated 1-100 keV luminosity. By normalising the 0.1-2 keV region of a power-law model 
with the same Ah to the entire HRI count rate of the source, we can derive this upper limit as 
L ~ 12 x 10 41 ergs _1 . The specific case of X-ray illumination of dense tori around AGN was 
modelled by Krolik & Lepp (1989), who found that ^i_os(i) ~ 10~ 4 Lx-ray for Seyfert Is and 
L 1 _qs(i) ~ 10~ 5 Lx-ray for Seyfert 2s. The 1-0 S(l) luminosity in a 5" aperture (similar to that 
used to estimate the X-ray flux from this source) around the SW nucleus in our line image is 
around 2 orders of magnitude larger than predicted for X-ray excitation alone. An additional 
argument against significant X-ray excitation is the offset of ~500 pc between the peak H2 and 
AGN positions. 

If the NE nucleus also hosted an AGN, similar arguments would apply. If, on the other hand, 
we assume the X-rays are associated with star- formation (eg supernovae or binaries), we should 
use the luminosities directly from Wang et al. (1997). The simplest models of H2 for this case 
are those of Lepp & McCray 1983. They predict that the total 1-0 S(l) luminosity is 0.25% of 
the X-ray luminosity, ie we would expect 0.6-3 x 10 5 L of 1-0 S(l) emission, whereas the total 
luminosity observed in 5" is 12.5 x 10 5 L Q . Thus X-ray excitation is expected to contribute 5-25% 
of the 1-0 S(l). 

4-2.3. Shock Excitation 

Shocks can be very effective at exciting H2 molecules and here we consider models of fast 
& slow J & C shocks from Burton et al (1990). These were originally tailored to the supernova 
remnant IC 433, but the general results have a much wider application. 

fast J shock (100-300 km s" 1 ) would dissociate the molecules, which then reform downstream 
on grain surfaces via formation pumping (Mouri & Taniguchi 1995). The observed spectrum 
would be similar to that of UV fluorescence as the molecules are created in energetic states 
which then decay radiatively (see Table 1 in that paper). These authors give the emissivity 
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(photons emitted per molecule formed) of 1-0 S(l) as ~0.02, requiring the shock(s) to pass 
through 1500 M yr" 1 near each nucleus if 50-70% of the fluorescent emission is due to 
fast J shocks. Additional tests are needed before we can confirm or rule out this process. 
That an equal flux of Br7 would be produced does not contradict the observations, as long 
as some other process can make up the bulk of the Br7 without producing significant H2 
emission. Observations of the J and H band [Fe 11] lines would be useful, as fast J shocks 
would destroy grains, enhancing the gas-phase Fe + abundance so that the line strength 
becomes comparable to 1-0 S(l) (Hollenbach & McKee 1989). 

slow J shock could be possible and would not dissociate the molecules. Shock velocities 
in the range v s = 5-15 kms" 1 result in strong H2 lines and a post-shock temperature 
T = (3900 K)(w s /10 kms" 1 ) (Shull & Draine 1987). However, in normal ionisation fractions, 
magnetic field strengths, and gas densities such shocks are expected to revert to C type 
(Hollenbach, Chernoff & McKee 1989). 

fast C shock with v s ~ 40 kms -1 heats the gas to > 2000 K (dependent on v s ) and produces 
strong H2 line emission, which could easily account for the bulk of the thermal excitation. 
We estimate a mass of 1600 and 1000 M Q of hot gas are needed near the NE and SW nuclei. 
The flow time through the hot shock structure of ~ 100 yr (Draine et al. 1983), then suggests 
that 10 and 16 M yr -1 of gas must be shocked to sustain the line's thermal luminosity. 

slow C shock is unlikely as the peak temperature is rather low, < 300 K, resulting in very weak 
H2 emission. This would contradict the strong observed emission and derived temperatures. 

We have considered the various mechanisms which might give rise to the H2 emission in 
the vicinity of each nucleus (R <^ lkpc), and find the simplest solution is that outlined below. 
The non-thermal emission is a combination of UV fluorescence (1-0S(1)/Br7 ~ 0.1) and fast 
100-300 kms -1 dissociative J shocks (1-0S(1)/Br7 ~ 1), the latter probably accounting for 50-70% 
with the ratios adopted here and heating ~1500MQyr _1 gas near each nucleus. The thermal 
emission is due in part to X-ray heating but mostly fast 30-50 kms -1 non-dissociative C shocks, 
heating only ~ 10-16 Mq yr -1 of dense clouds. In terms of the 1-0 S(l) line, the contributions are 
approximately: for the NE nucleus 5-25% X-ray heating, <5% UV fluorescence, 15-20% fast J 
shocks, 55-75% fast C shocks; for the SW nucleus <20% UV fluorescence, 10-30% fast J shocks, 
70% fast C shocks. 

4.3. Spatial Analysis 

In order to reach sufficient signal-to- noise for the above analysis, we examined only the 
integrated flux around the nuclei and ignored the morphology. A spatial anaylsis of the brighter 
features along the slit can also be fruitful, as is apparent from Fig |6|. The upper row shows the 
continuum, 1-0 S(l) line and Br7 line fluxes along the slit from the SW nucleus (left) to the NE 
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nucleus (right). That the (relative) line fluxes do not exactly tally with those given in Table |l| is 
due to different apertures: the data in the table were extracted from a 5" circular aperture, those 
in the figure from 0.61" x 1.22" boxes. Also, no structures less than 1-1.5" will be visible, so the 
details seen in the images in Fig |l] are not apparent. 

The lower row of the figure is more interesting: it shows that the ratios of 1-0 S(l) to both 
Br7 and the continuum increase away from the nuclei. Thus although the H2 emission right in 
the nuclei may be associated with Br7 and the 2/mi continuum (eg via star formation processes), 
that more than ~1" away is not. Remarkably the ratios remain high even 2-3" away, a physcial 
distance of 1-1.5 kpc, although the surface brightness is very low. 

4.4. Extinction 

The extinctions derived from the Ha/H/3 ratio in Veilleux et al. (1995), and also using the 
Ha/Br7 ratio, are very small (Ay^3, Table [l]). In the case of the NE nucleus, they are almost the 
same and suggest a correct measure for the total extinction. However, for the SW nucleus, the Ay 
found using the Br7 line is significantly larger, an indication that we are probing to greater optical 
depths with the less-absorbed near-infrared line. Additionally continuum fits to the spectra give 
Ay ~ 4, remarkably consistent (expecially in the K-band where Ak = O.lAy) given the way they 
were determined, and that they probe the stellar continuum rather than the excited gas. Lastly, 
we measured the Paa/Br7 ratio, which yielded somewhat higher values of Ay ~ 7 for the NE 
and Ay ~ 13 for the SW nuclei; we must disregard these since the Paa lines are seen at very low 
(<40%) atmospheric transmission, and systematic errors of only 10-30% would bring them in line 
with the other estimates. 

5. Discussion 

In this section we draw together the various lines of evidence that we have examined, in an 
attempt to paint a consistent picture of the events occuring in Mkn266. 

5.1. North East Nucleus 

The NE nucleus is classified from its optical spectrum as a LINER by Osterbrock & Dahari 
(1983) and as a Seyfert 2 by Veilleux et al. (1995). Some confusion is unavoidable as explanations 
for LINER spectra are typically either weak AGN or shock excitation, maybe due to a starburst 
(Filippenko 1996). Both we (at 2/im) and MGAH (at 2-20 cm) have found that this nucleus 
is resolved on scales of several hundred parsecs, making an AGN an unlikely proposition. 
Additionally, Smith et al. (1998) found evidence for simple structure at VLBI 18 cm scales. The 
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flux they detected was 4.5 mJy, only 20% of that in the entire nucleus, suggesting that even if 
there is an AGN component, it certainly cannot dominate the power output. Instead we assume 
that all the Br7 is due to recent star formation, and using models of Leitherer et al. (1995) we 
estimate the supernova rate (z^sn), K-band luminosity (-&k)> and bolometric luminosity (Lbol) with 
which it is likely to be associated. We employ models with solar metallicity and a Salpeter (slope 
2.35) initial mass function (IMF) in the range 1-100 M . 

If active star formation occured over a short period and has since evolved passively, then 
it must be very young as Br7 flux and equivalent width fall off very quickly with time; the 
directly measured Br7 equivalent width (Wbt-j) of 10 A (in a 2" aperture) provides an upper limit 
to the age. This is an unrealistic number to use since at least some of the continuum is from 
an old underlying population (as evidenced by the extended continuum): decomposition of the 
deconvolved image suggests that in the central 2" , roughly 70% of the continuum is from the core 



and the rest extended. Additionally in Section iA_ we showed that supergiants make up ~25% of 
the continuum. Taking these two cases as the extremes implies that the Wsr-y associated with 
the most recent episode of star formation should lie in the range 15-40A. This constraint give 
a range of 4.5-6.5 Myr, consistent with the CO index CO p h = 0.12. Such a young age implies 
3-8x10 7 Mq of stars would have had to form effectiv ely instantaneously (<^ 1 Myr) in order to 
reproduce the observed Br7 flux. We would then predict the parameters z^sn = 0.03-0.08 yr _1 , 
L K = 2-7 x 10 8 L , and L hol = 3-6 x 10 10 L Q . 

Models for continuous star formation may also be appropriate for this object, and a plausible 
scenario might invoke 'punctuated star formation' in which clusters of 10 5 -10 6 M formed at 
different times over this period. Wbt^ = 15-40 A then implies an age of 60-500 Myr (or more), 
during which entire period the stars have formed with an average rate of 1.7M yr~ 1 (ie, a few 
clusters every Myr). 60 Myr is the time at which the maximum supernova rate i^sn = 0.03 yr _1 
is reached and thereafter sustained. Unfortunately the CO index does not constrain the age 
range, although both Lk and Lboi should since they increase slowly with time as more late-type 
giants are accumulated. Yet even after 500 Myr the model predicts only Lk = 5 x 10 s L and 
L bol = 4x 1O 1O L . 

The parameters for these two models are summarised in Table ||| as Models l and 2 
respectively, and are compared with observed quantities. For these we used the K-band luminosity 
from this paper and, as above, suppose that somewhere in the range 25-70% of that in 2" comes 
from the starburst. The bolometric luminosity is the infrared (8-1000 /jm) luminosity taken from 
Sanders et al. (1991). We assume that each nucleus contributes in rough proportion to its Br7 
and X-ray fluxes, amounting to 30-50% for the NE nucleus. The model predictions for these two 
quantities are factors of 3-4 less than those observed. 

The supernova rate ^sn has an even larger discrepancy. To estimate it from the 20 cm 
continuum we used the relation derived by Condon & Yin (1990), between non-thermal radio 
continuum and supernova rate. By considering independent estimates of these for the Galaxy they 
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found: 

LntIWHz- 1 ] ~ 1.3 x 10 23 (i/[GHz])-°- 8 ^[vr" 1 ] 

This is very similar to the relation derived for M82 by Huang et al. (1994), who find a coefficient 
of 1.1 ± 0.5 x 10 23 for that galaxy by considering the cumulative number of sources as a function 
of increasing diameter. Additionally, by using the temporal behaviour of the bright SN1979c 
in M 100 as a template, Colina & Perez-Olea (1992) derived a coefficient of 0.9 x 10 23 for the 
same equation. So this relation appears to be quite robust, whether applied to the Galaxy or to 
the rather different environment of starbursts. From the 20 cm radio continuum (MGAH) which 
has a spectral index of 0.6 and hence will be dominated by the non-thermal component at this 
wavelength, we derive i^sn = 0.45 yr _1 . 

These calculations provide a strong indication that the starburst scaling obtained from the 
Br7 flux is too small by an order of magnitude. There are potentially several ways around this 
problem. A number of authors have suggested that the IMF might be truncated or have a steep 
slope, both of which reduce the number of most massive stars. These two cases were also modelled 
by Leitherer et al. (1995), and W Brj = 20 A gives maximum ages of 30Myr (Miller-Scalo IMF 
with slope 3.3) and HMyr (upper mass limit set to 30 M Q ). Although the normalisation of the 
model to the observed ionising flux then increases, the predicted supernova rate actually decreases. 
For the truncated IMF it is because the age is so small that many of the progenitor stars are too 
young to produce supernovae; for the steep IMF it is mainly because there are fewer supernova 
progenitors. 

The required effect of reducing the observed Br7 flux could also be achieved if there was dust 
internal to the ionised nebulae, which then competes for Lyman continuum (Lyc) photons. Perhaps 
more likely is that if the molecular clouds have been shredded by multiple supernovae, the nebulae 
may be density bounded so that many of the Lyc photons escape. With the measured Ay ~ 2 
more than 90% of these photons would be mopped up by dust in the ISM, contributing to the 
far-infrared flux without increasing the observed UV. In either case the hydrogen recombination 
line strengths would be reduced by the same factor. An intriguing side-effect is to alter the 
ratio of Ha or H/3 to the high ionisation species such as [Om] . Since lines with higher ionisation 
potentials are produced deeper in the nebula, their strengths will not be affected if the nebula 
is density bounded (except in rather extreme cases). Such an effect could account for the high 
ratio of [Om]/H/3 measured by Veilleux et al. (1995) leading to their classification as a Seyfert. 
We can hypothesise a 'intrinsic' Br7 flux (ie the one we ought to measure if all the Lyc photons 
ionised Hi), which would be many times larger than that observed. However, this does not solve 
our problem since Wbt^ must also increase by the same factor, forcing us again to adopt a very 
young starburst age. 

One remaining hypothesis involving star formation is that at some point 10-30 Myr ago 
a catastrophic event may have triggered a single highly intense burst of activity, forming 
~ 7 x 10 8 M during this period. The starburst would now have evolved enough that it contributes 
little to the current Br7 luminosity (1-2%) while accounting for the radio continuum, Lk, and 
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Lboi- This is shown as Model 3 in Table |3[ The corollary is that this leaves almost no age 
constraints on the star formation responsible for the Br7. However, the difficulty here is that the 
radio and 2 /im continua trace out different regions so it seems unlikely they could originate from 
the same episode of star formation. 

A different approach is to consider an alternative origin for the synchrotron radiation - which 
must nevertheless arise from electrons accelerated by shocks since we have shown that an AGN 
is not a proposition. Perhaps the two most important pieces of evidence are: (1) the 1-0 S(l) 
and radio continuum appear to be related to each other because of their shock origin, and (2) 
although the radio continuum is resolved, both it and the 1-0 S(l) are very compact (on scales 
of 0.4") - and, crucially, more so than the 2 (im continuum (0.7"). Our spectra show that the 
2 /im continuum is dominated by late-type stars, and presumably traces out the region of recent 
star formation. We argue that perhaps the radio and 1-0 S(l) are not directly associated with 
the star formation. Although the spatial scales are too small to make a definitive statement, the 
morphology does appear to be similar to that of NGC 6240 where the gas has settled between 
the two continuum nuclei and is radiating strong thermal 1-0 S(l) emission (van der Werf et al. 
1993, Sugai et al. 1997b, Tacconi et al. 1999) . If NGC 6240 were moved away so that its nuclei, 
instead of being 1.5" apart, were separated by only 0.5", it would begin to look like Mkn266 NE. 
We speculate whether Mkn266 NE could itself consist of two or more (unresolved) continuum 
nuclei with gas settling between them. As the models we have considered imply, star formation 
has occured recently in these nuclei, but the supernova rate is very low so essentially all the radio 
continuum originates in the settling gas. A serious difficulty with such an interpretation is that 
Mkn266 NE has a 3.6cm intensity (scaled from the 2 cm, MGAH) of 7.0 mJy within 0.3"x0.4", 
while NGC 6240 has diffuse 3.6 cm emission at the level of 0.2 mJy in a 0.23" x 0.32" beam (Colbert 
et al. 1994) - some 35 times fainter. Since we cannot invoke an AGN, it is hard to conceive of a 
method to produce so much more non-thermal radio emission without reverting to fast shocks. 
Bell (1978) showed that the radio emissivity of shocked gas is strongly dependent on the shock 
velocity v s : 
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So for a typical spectral index of a = 0.8 we find that e(u) oc v}.^ for v s ^ 2000 km s^ 1 ; while 
e{u) oc vf 2 for v s ;> 7000kms _1 . To produce such intense synchrotron emission requires shocks 
with speeds of several thousand kms -1 , 3-10 times faster than the 300kms _1 seen in the CO line 
in NGC 6240 by Tacconi et al. (1999). However, in NGC 6240 higher velocities than in the CO line 
were observed in 1-0 S(l) which has a FWHM of 550 km s" 1 and FWZI of 1600 km s -1 (van der 



Werf et al. 1993). Also, in Section 4^ we saw that the fraction of H2 excited by dissociation due to 
fast shocks is much higher in Mkn266 than in NGC 6240; and to account for this ~100 times more 
gas in Mkn 266 is shocked at high velocities than at the low velocities expected in dense molecular 
clouds. Bell's equation shows that, keeping all other parameters the same, increasing the shock 
velocity from 300kms _1 to 3000 km s -1 would increase the synchrotron intensity by a factor of 



-16- 



30, similar to the observed difference between NGC6240 and Mkn266. It does not seem totally 
unreasonable then to extrapolate from NGC 6240 to a more extreme environment in Mkn 266 NE: 
if gas has settled into a disk structure on scales several times smaller than in NGC 6240, the 
velocities required to support it centrifugally will be many times greater and this could give rise to 
the fast shocks required to produce the intense synchrotron emission. High resolution observations 
of the CO in this nucleus will be needed to decide whether this speculation is indeed borne out. 

One final point to consider is the existence of faint but very extended wings out to at least 
1-1.5 kpc in the 1-0 S(l), as indicated in Fig ^. The lower limits to the S(l)/Br7 ratio and the 
high S(l) equivalent width indicated in Fig ^ imply an origin not associated with star formation. 
A plausible alternative is that the lines are excited by shocks which are driven into clouds in the 
turbulent medium of the X-ray halo. 

5.2. South West Nucleus 

This nucleus has also been identified as a LINER (Veilleux et al. 1995) and Seyfert 2 
(Osterbrock &; Dahari 1983). The classification as a LINER rested solely on the ratio 
[Olll]/H/3 = 1.4 because the diagrams used to distinguish LINERs, Seyferts, and starbursts all 
plotted this value, against the ratios of [Nil], [Sn], and [Oi] to Ha. Evidence for an AGN comes 
from MGAH who observed an unresolved radio core in the SW component using a beamwidth of 
0.3" x 0.4", and our deconvolved continuum image. An additional line of reasoning put forward 
by Wang et al. (1997) is that photoionisation by the hard radiation field of an AGN is needed 
to explain both the luminosity of the large-scale optical emission line nebula as well as the line 
ratios. As discussed in Section |H based on the 2 /jm continuum we would argue for a combination 
of Seyfert plus starburst. 

The surrounding emission, extending linearly about 1" both north and south (radio as well 
as infrared) would be due to star formation. The Bi~7 follows the distribution to some extent, 
although the resolution is rather poor; for example, the extension to the east is apparent in both 
continuum and line images and may be due to an extranuclear Hn region. Unfortunately, a 
quantitative analysis is not possible with the current spatial resolution, since we cannot properly 
separate the AGN and starburst contributions. 

More interesting is the H2, which peaks to the north east between the galaxies in a region 
without any associated Br7 or continuum, and 500 pc from the nucleus. This would seem to rule 
out any connection with either the AGN or star formation. In Section we showed that 70% 
of the 1-0 S(l) emission in this nucleus was thermally excited and deduced that it arose mostly 
in ~ 40kms~ 1 C shocks. Together these suggest that it may due to shocks resulting from the 
interaction being driven into clouds. Both we and Veilleux et al. (1995) found a radial velocity 
difference between the nuclei of 50kms _1 , setting a lower limit on the speed of the interaction. It 
is likely that the actual speed is faster, but shock speeds similar to that required can be produced 
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when a high speed shock in a low density medium encounters a high density medium, the process 
favoured for the H 2 excitation in NGC 6240 by van der Werf et al. (1993) and Sugai et al. (1997b). 
However, it is not possible to comment here on whether it is due to the global interaction, or if it 
results from the gas dynamics more closely associated with the SW nucleus. 

5.3. The Central Radio Continuum 

One of the puzzling observations of Mkn 266 is the non-thermal radio continuum midway 
between the nuclei, stronger than the extended emission in the SW nucleus and accounting for 
more than 1/4 of the total 20 cm emission (MGAH). What makes this so surprising is that there 
is no evidence for any other enhanced emission of any sort at this position. 

MGAH argued that it could result from shocks as clouds collide near the interface of the gas 
disks. Our most severe observational constraint on this model is the non-detection of 1-0 S(l) 
in the locale, at a 5a limit of 10~ 18 Wm~ 2 (in a 5" aperture, similar in size to the 3"x7" radio 
emission). The MGAH scenario is ideally suited to the interacting galaxy model of Jog and 
Solomon (1992), which describes pairs of gas rich spiral galaxies colliding at 300 km s -1 . This 
model has an important addition to previous models (eg Harwit et al. 1987) which have invoked 
collisions between GMCs to explain the observations. The authors note that GMCs have a low 
volume filling factor /gmc = 0-01 and so are unlikely to collide; it is the Hi clouds with /m = 0.1 
which collide. The GMCs are assumed to have a core radius of 25 pc with an average density 
jih 2 = 100 cm~ 3 ; the Hi clouds are given a radius of 5pc and an average density wh = 20cm -3 , 
but with the larger filling factor, the total mass is the same as the H2. 

The fast shocks in the Hi clouds are slowed in the denser GMCs to ~ 50kms -1 , which are 
still nevertheless highly supersonic; the compression timescale of ~ 2 x 10 4 yr then implies that in 
total only about 10% of the GMC gas is shocked. The cooling timescale is on the order of 100 yrs 
(Drain, Roberge & Dalgarno 1983), so only 0.5% of the shocked H2 is radiating at any time. 

At a distance of 115 Mpc, and taking the line-of-sight dimension as an average of the others, 
the volume over which the radio continuum is emitted is 24kpc 3 . The H2 mass inside this would 
then be 10 9 M , of which 5 x 10 5 M would be hot - roughly the same amount as in NGC 6240 
(Sugai et al. 1997b), and about 500 times more than is implied by our 5a upper limit to the 
1-0 S(l) line. Even allowing for large uncertainties in the parameters and timescales, invoking 
shock excitation in clouds induced by the interaction appears unlikely. 

We offer a slightly different interpretation. Wang et al. (1997) have argued that X-ray and 
optical emission line observations point towards a radial outflow of gas at several hundred kms -1 
which has swept up much of the ISM, and they conclude that the NE starburst and SW AGN 
play comparable roles in photoionising the nebula. At distances of a few kpc, the outflow velocity 
is ~ 300 km s -1 . Our hypothesis is that the radio continuum traces the interface between the 
expanding shock fronts from the two nuclei. When added to the relative velocity of the galaxies, 
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the relative velocities of these shocks will be approaching 1000 km s . Copious quantities of 
non-thermal radio emission will be produced even if the density of the shocked material is fairly 
low, but there will be no 1-0 S(l) and little hydrogen recombination flux as the gas is heated 
to 10 6 -10 7 K. A test of this would be to observe the region with high spatial resolution in the 
0.1-2 keV range, to see if there is a local increase in the temperature of the X-ray emitting gas 
over the ambient temperature in the bubble of 10 6 K. 

6. H2 Excitation: The Wider Perspective 

The classical case of H2 emission occuring between the nuclei in a merger is that of NGC 6240 
(Herbst et al. 1990, van der Werf et al 1993). The H2 level population is consistent with solely 
thermal excitation and, combined with the absence of [Fen] in the same region, indicates that the 
observed emission arises from non-dissociative shocks, occuring as a result of the merging process 
(Sugai et al. 1997b). Recent CO radio observations (Tacconi et al. 1999) have confirmed the 
presence of 2-4 x 10 9 M of molecular gas concentrated here. 

Other radio CO observations of a number of galaxies at 1-3" resolution by Gao et al. (1997) 
tend to support this more generally, and Scoville et al. (1997) argue that 2/3 of the molecular gas 
in Arp 220 relaxing into such a disk is consistent with scenarios in which gas in merging systems 
settles into the centre faster than the stellar nuclei. They argue that efficient removal of angular 
momentum from the observed dense centrifugally supported disk could play an important role in 
the evolution of powerful starbursts/AGN. They estimated a mass of 5.4 x 10 9 Mq in the disk 
with a mean density of 2 x 10 cm -3 . Based on the model of Jog & Solomon (1992) discussed 
in Section |5,3| , we estimated a ratio of cold to hot (shocked) molecules of 2000. This is not so 
different from the ratio of found in NGC 6240: Tacconi et al. (1999) estimated the cold H2 mass to 
be 2-4 x 10 9 M , while Sugai et al. (1997b) found a cold H 2 mass of 2 x 10 5 M , yielding a ratio of 
10 4 . Even for this case the 1-0 S(l) emission should be easily observable. For example, for a merger 
at 200 Mpc to produce a 1-0S(1) flux of 10" 18 Wm~ 2 , an H 2 disk mass of only ~ 2 x 10 7 M is 
needed. Given the examples of Arp 220 and NGC 6240, this seems a perfectly feasible proposition 
for any of the luminous infrared galaxies (almost all at distances < 200 Mpc) observed by Sanders 
et al. (1991), in which the CO luminosity implies total H2 masses in the range 1O 9 -1O 1O M . 

We have now observed several other such systems similar to NGC 6240 by targeting close 
mergers, although with current observations it is not possible to say whether the gas has formed a 
disk structure in these galaxies. These include Mkn551 (Sugai et al 1997a), Mkn266 (here), and 
Mkn496 (Sugai et al. 1999a). Also, the line images of NGC 3256 (Kotilainen et al 1996) suggest 
that there is an enhancement of 1-0 S(l) between the two nuclei with respect to both the Br7 and 
[Fen], which fall to a minimum. This phenomenon is not limited to the luminous infrared galaxies, 
and has also been observed in some blue compact dwarf galaxies (II Zw 40, NGC 5253, and He 2-10, 
Davies et al. (1998)) where again the 1-0 S(l) bears little or no resemblence to the Br7 or 2fim 
continuum. However, the 1-0 S(l) surface brightness is low in these interacting systems, and 
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it could only be detected because they are all nearby (<10Mpc). The more luminous galaxies 
with correspondingly larger gas masses do form a better sample for further study, and although 
such emission is not always observed (eg Arp299, Sugai et al 1999b), it may nevertheless be a 
relatively common phenomenon and should be observable with the current generation of infrared 
instrumentation. 

7. Summary & Conclusions 

We have presented near infrared data on Mkn 266 consisting of 2/xm continuum, and Br7 and 
1-0 S(l) emission line images, as well as K-band spectra. These are analysed in conjunction with 
data from the literature, primarily the radio continuum (MGAH). 

There are 3 main observations to note from the images. 

(1) In the NE (LINER) nucleus the 1-0 S(l) is similar in size (^0.45") to the radio continuum, but 
more compact than the 2/txm continuum (0.7"). 

(2) In the SE (AGN+starburst) nucleus, the 1-0 S(l) is offset from the continuum emission by 0.9" 
(500 pc), and has no associated Br7 emission or radio continuum. 

(3) Midway between the two components where there is a region of strong radio continuum, we 
find no enhanced near infrared line or continuum emission at all. 

From our analysis of the spectra we note the following. 

(1) By fitting stellar templates to parts of the spectra away from emission lines, we are able to 
remove all the absorption features, and hence measure the fluxes of even weak emission lines. 
Simultaneously we find that 20-25% (NE) and 50% (SW) of the continuum is due to supergiants, 
the rest originating in late type giant stars. 

(2) Level population diagrams of the hot H2 molecules indicate that in the NE and SW nuclei 
respectively, 81±2% and 71±12% of the 1-0 S(l) is thermally excited (to temperatures of 1500±60 
and 2460±410 K), mostly by ~40kms _1 C-shocks. Additionally, much of the non-thermal emission 
may arise in faster dissociative J-shocks. Pure thermal models (with either one or two components 
of different temperatures) are effectively ruled out. 

(3) Although the 1-0 S(l) surface brightness is very low away from the nuclei, its ratio to Br7 and 
the continuum remains high even out to 2-3". 

The conclusions we draw are listed below. 

NE nucleus: the extended nature argues against an AGN, but nor can any combination of 
star formation scenarios simultaneously account for the observations. Both the morphologies of 
the 1-0 S(l) and 2/jm continuum, and the (predominantly) shock origin of the H2 emission and 
radio continuum, are reminiscent of NGC 6240. We speculate that this nucleus may resemble 
NGC 6240 but on a smaller physical scale, and with the higher velocities that would necessarily 
occur in a smaller centrifugally supported disk. 
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SW nucleus: the morphologies suggest an AGN with circumnuclear star formation, but a 
quantitative analysis is not yet possible due to the small angular scales involved. Because the 
offset 1-0 S(l) has no other associated emission it cannot originate in star formation, and may 
instead be due to shocks driven into clouds as a result of the interaction. 

intermediate region: The lack of any enhanced 1-0 S(l) emission suggests that the radio 
continuum cannot trace out regions shocked by the interaction. We propose that it indicates the 
interface between the expanding winds from each nucleus, where combined shock velocities would 
be high enough (~1000kms _1 ) that neither 1-0 S(l) nor hydrogen recombination lines would be 
observed. 
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Tablc 1. Mkn266: Basic Quantities 



Quantity 


Unit 


North-East Nucleus 


South- West Nucleus 






5" 


2" 


5" 


2" 


2 /im continuum a 


[mJy] 
[mag] 


9.43 ± 0.06 
12.10 ±0.01 


4.09 ± 0.03 
13.00 ±0.01 


14.07 ±0.06 
11.67 ±0.01 


5.76 ± 0.03 
12.63 ±0.01 


Br7 line a 


[10- 18 Wm- 2 ] 


3.4 ±0.4 


2.6 ±0.2 


5.7 ±0.4 


2.8 ±0.2 


1-0 S(l) line a 


[10- 18 Wm- 2 ] 


3.2 ±0.2 


2.2 ±0.1 


2.3 ±0.2 


0.7 ±0.1 


1-0S(1)/Br 7 a 




0.94 ±0.13 


0.85 ± 0.08 


0.40 ± 0.05 


0.25 ±0.04 




[A] 


5.7 ±0.7 


9.9 ±0.8 


6.3 ±0.4 


7.6 ±0.5 


Av b 


[mag] 




2.3 


0.9 




A V C 


[mag] 


2.7 ±0.2 


2.1 ±0.2 



a Measured in aperture given; errors are statistical, determined from noise in the images 
b from Ha/H/3 ratio, Veilleux et al. (1995) 

c from Ha/Br7 ratio, where Ha is from Veilleux et al. (1995) and Br 7 is measured in same region 
(3.7" length section of a 2" wide slit at position angle 36°). Errors are from Br 7 flux only. 
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Table 2. relative line strengths & derived quantities' 



Line 


Wavelength 


E u /k h 


A„i b 

x Ul 


Relative strength 






K 


10- 7 s -1 


NE C 


SW d 


1-0S(3) C 


1.9576 


8365 


4.21 


0.76 




1-0 S(2) 


2.0338 


7584 


3.98 


0.33 


0.46 


2-1 b{6) 


2.0735 


13890 


5.77 


0.12 


0.22 


1-0 S(l) 


2.1218 


6956 


3.47 


1.00 


1.00 


2-1 S(2) 


2.1542 


13150 


5.60 


U.U ( 




3-2 S (3) 


2.2014 


19086 


5.63 


0.04 




1-0 S(0) 


2.2235 


6471 


2.53 


0.27 


0.34 


2-1 S(l) 


2.2477 


12550 


4.98 


0.12 


0.24 


1-0 Q(l) 


2.4066 


6149 


4.29 


0.81 


0.65 


1-0 Q(2) 


2.4134 


6471 


3.03 


0.28 


0.35 


1-0 Q(3) 


2.4237 


6956 


2.78 


0.69 


0.85 


Paa f 


1.8756 






5.32 


15.00 


Br7 


2.1661 






0.54 


1.82 


He I 


2.0587 






0.34 


0.85 


a Errors are 


estimated from the 


residual continuum 


near each line. 


Due to poor 


atmospheric 


transmission and continuum subtraction shortwards of 1.9/im, additional errors are present in this 


region and we 


are unable to measure the Br<5 (1.95/im) 


or 1-0 S (4) (1.89//m) lines. 




b For H 2 


lines only; 


taken from 


the UKIRT 


Astro-Utilities page: 



http:// www.j ach.hawaii.edu / J ACpublic /UKIRT /astronomy /| 



c Ratios are normalised to 1-0 S(l); la errors are ~0.02 

d Ratios are normalised to 1-0 S(l); lex errors are 0.05-0.07, except 0.11 for the Q-branch 

e As there is no evidence for an AGN in the continuum, we have assumed there is negligeable 
[Sivi] emission at 1.96/um. 

f As explained in a, there may be a large (10-30%) systematic error on this line making extinction 
estimates highly uncertain. 
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Table 3. Mkn266 NE: Star formation models a 



Quantity 


Unit 


Model l b 


Model 2 b 


Model 3 C 


Observed d 


Note d 


Age 


10 6 yr 


4.5-6.5 


60-500 


20 






Stellar mass 


10 7 M 


3-8 


10-85 


70 








yr" 1 


0.03-0.08 


0.03 


0.40 


0.45 


from 20 cm radio 


L K 


1O 8 L 


2.1-6.9 


2.7-6.0 


20.6 


6-18 


25-70% of flux in 2" 


Lbo\ 


1O 1O L 


3.3-5.7 


2.9-3.8 


12.0 


9-15 


30-50% of total L m 



a Using the models of Leitherer et al. (1995) 

b Instantaneous (1) and continuous (2) models normalised to Br7, assuming that 25% (lower ages) 
and 70% (higher ages) of the continuum within 2" is from the starburst. 

instantaneous model (3) scaled to i^sn, Lk, and Lboi- 
d The quantities measured, and where they came from 
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1.7 kpc 




17 kpc resolution 




J I I I I I I I I I I I I I I I I I L 



1.7 kpc resolution 




Fig. 1. — Contour maps of the 20" region around Mkn266. Upper left: 2 /jm continuum. Contours 
are from 2<t pixel -1 at intervals of 6a pixel -1 , where la pixel -1 = 0.072 mJyarcsec -2 . Upper 
right: 2 /im continuum after rebinning and deconvolving with 30 iterations of the Lucy algorithm 
(see Section 1 for details) at arbitrary linear intervals. Lower left: 1-0 S(l) line superimposed 
on continuum greyscale. Contours are from 2.5crpixel -1 at intervals of 2. 5a pixel -1 , where 
la pixel -1 = 1.6 x 10 -19 W m 2 arcsec 2 . Lower right: Br7 line (smoothed using a 3 x 3 median 
filter, with little effect on the resolution) superimposed on continuum greyscale. Contours are from 
2.5<rpixel -1 at intervals of 1.5a pixel -1 , where lapixel -1 = 2.1 x 10 -19 Wm -2 arcsec -2 . 
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Fig. 2. — Left: Measured profiles of a star, 2 fim continuum and 1-0 S(l) line in the North-East 
nucleus. Azimuthally averaged intensity profiles are derived from flux in concentric annuli. Right: 
ratio of 1-0 S(l)/2 fxm. (ie normalised equivalent width of the line), decreases as a function of radius. 
These diagrams show clearly that the core of the 1-0 S(l) is (a) unresolved, whereas the continuum 
is marginally resolved; and (b) has faint wings of emission at 1-2". 
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Fig. 3. — Integrated spectra of the NE (upper) and SW (lower) nuclei of Mkn 266. The important 
emission and absorption features are indicated. Fitting and subtraction of the continuum with 
stellar teplates is important in order to account for absorption features close to emission lines (eg 
the H2 Q-branch). 
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Fig. 4. — Left: Population diagram for the hot H2 molecules in the NE nucleus, derived for 11 
line transitions. The crosses denote the combined fit to the data of a thermal component (dotted 
line, with derived temperature T = 1500 K) and a non-thermal component (Black & van Dishoeck 
(1987) UV fluorescent model 14). The reduced chi-square Xu = 0.89, close to the expectation value 
(0.92) for 8 degrees of freedom. Right: Confidence regions for the thermal excitation temperature, 
T, and the non-thermal fraction of 1-0 S(l), fuv- The la uncertainties give T = 1500 db 60 K and 
fuv = 0.19 ± 0.02 (ie fluorescent excitation accounts for 57% of the total H2 line emission). 
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Fig. 5. — Left: Population diagram for the hot H2 molecules in the SW nucleus, derived for 8 
line transitions. The crosses denote the combined fit to the data of a thermal component (dotted 
line, with derived temperature T = 2460 K) and a non-thermal component (Black & van Dishoeck 
(1987) UV fluorescent model 14). The reduced chi-square xt = 1-08, close to the expectation value 
(0.87) for 5 degrees of freedom. Right: Confidence regions for the thermal excitation temperature, 
T, and the non-thermal fraction of 1-0 S(l), fuv- The la uncertainties give T = 2460 ± 410 K and 
fuv = 0-29 ± 0.12 (ie. fluorescent excitation accounts for 68% of the total H2 line emission). 
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Fig. 6. — Upper: 2/xm continuum, 1-0 S(l) flux, and Br7 flux, traced for 40" along the slit across the 
SW and NE nuclei. In the latter 2 figures, the vertical bar indicates teh SW continuum centroid; 
unlike the images, no offset is seen between the 1-0 S(l) and 2/xm continuum due to the poorer 
spatial smapling and seeing. Lower: ratio of 1-0S(1)/Br7 and l-0S(l)/2/im (equivalent width). 
Away from the nuclei, both ratios increase. In the NE nucleus, the initial decrease is the same as 
that seen in Figure 2. 



